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ABSTRACT 

This  thesis  presents  the  application  of  EIBENVALUE 
SENSITIVITY  ANALYSIS  and  SINGULAR  VALUE  ANALYSIS  to  the 
control  o+  a  SEA-SKIMMING  supersonic  missile, in  the  vertical 
plane. 

The  study  is  divided  in  -Four  basic  parts: 

a)  The  development  o-f  the  model. 

b)  EIGENVALUE  SENSITIVITY  ANALYSIS  with  respect  to  the 
variation  o-f  the  aerodynamic  parameters  o-f  the 
autopilot/air-frame  o-f  the  missile. 

c)  Analysis  o-f  the  time  response  with  respect  tc  the 
variation  o-f  the  aerodynamic  parameters. 

d)  Robustness  analysis  and  improvement  o-f  the  system, 
using  the  SINGULAR  VALUE  ANALYSIS. 

All   the   analysis   is   based   in   results  o-f  simulation 

programs    using    the    so-ftware   available   at   the   NAVAL 

POSTGRADUATE  SCHOOL. 
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I.  INTRODUCTION 

The  majority  o-f  tactical  missiles  used  against  surface 
targets  are  SEA-SKIMMIN6  and  -fly  as  close  as  possible  to  the 
surface  of  the  water  in  order  to  make  difficult  the 
detection  and  reaction  of  the  enemy. 

This  work  addresses  the  problem  of  analysing  the  design 
of  the  altitude  control  system  of  a  supersonic  SEA-SKIMMING 
missile  in  both  aspects  of  SENSITIVITY  and  ROBUSTNESS,  with 
respect  to  the  variations  of  the  aerodynamics  parameters. 

Those  variations  can  develop  from  uncertainties  in  the 
model  as  well  as  from  changes  in  the  flight  conditions  of 
the  vehicle. 

The  SENSITIVITY  analysis  makes  it  possible  to  verify, 
for  a  given  model  which  parameters  are  more  important  for 
the  desired  flight  path. 

The  ROBUSTNESS  analysis  has  the  objective  of  verifying 
if  the  system  could  be  affected  by  perturbations  or  noise 
and  will  be  the  base  on  calculations  of  modified  feedback 
gains  which  are  capable  of  producing  the  desired  flight  of 
the  missile  despite  perturbations. 

The  contents  of  the  different  sections  are  as  following: 

Section  II  presents  the  model  and  the  state  equations  of 
the  system  as  well  as  the  calculations  to  complete  the  chain 
of   control   of  the  missile,  considering  the  restrictions  of 
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the  SEA-SKIMMING  scenario,  where  the  altitude  control  system 
cannot  tolerate  any  overshoot  to  avoid  that  the  missile  hits 
the  water. 

Section  III  develops  the  SENSITIVITY  ANALYSIS  taking 
into  consideration  the  variation  of  the  eigenvalues  and  time 

response  o-f  the  system  o-f  state  equations  -for  variations  o-f 

the   aerodynamic   parameters   in   a   range   o-f   +25%  of  the 

original  value. 

In  section  IV  is  described  the  ROBUSTNESS  ANALYSIS  based 
upon  the  minimum  singular  value  of  the  return  difference 
matrix  as  function  of  the  frequency;  section  V  introduces  a 
design  technique  using  optimization  routine  in  order  to 
obtain  a  "complete  robust"  solution  and  section  VI  presents 
the  final  conclusions. 
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II.  MATHEMATICAL  MODEL  OF  THE  MISSILE 

The  work  assumes  a  SEA-SKIMMING  missile  with  the  -flight 
in  the  vertical  plane  according  to  the  path  shown  in  Figure 
II.  1. 
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Figure  II. 1  -  Missile  Trajectory  in  the  Vertical  plane 

The   missile   is   assumed  to  -fly  at  constant  velocity, at 
150  ft,  when  it  is  commanded  to  a  final  altitude  of  30  ft. 
Only  the  pitch  channel  will  be  analysed. 
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A.   THE  MISSILE 

The  missile  will  be  considered  as  presented  by  Arrow 
CRe-f.13.  This  model  has  been  selected  in  order  to  avoid 
classi-f  ication  problems.  The  geometry  as  well  as  its 
characteristics  Are    presented  in  Appendix  A. 

The  states  that  were  considered  in  our  model  are  de-fined 
in  Table  I  as  well  as  the  relevant  aerodynamics  parameters. 

In  the  Figure  1 1.2  we  have  the  schematic  representation 
o-f  the  uncoupled  pitch  control  chain  of  the  missile, 
neglecting  the  yaw  and  roll  movements; 


H 


Figure  II. 2  -  Block  Diagram  of  the  Pitch  Control   Chain    of 
the  Missile. 


where 

^zcb 

^c 


-  "bias"  acceleration; 


-  commanded  acceleration; 


achieved  acceleration; 
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v(t)  -  velocity; 
H<t)  -  altitude; 

H     —  commanded  altitude. 

s 

The  Transfer  Function  G<s)  represents  the  missile  itsel-f 
and  the  other  blocks  are  the  controller.  The  input  to  this 
block  will  be  the  "commanded"  acceleration  and  the  output 
Mill  be  the  achieved  acceleration. 

The  commanded  acceleration  comes  from  the  -feedback  of 
the  states  of  the  controller  (compensated  state,  velocity 
and  altitude)  added  to  a  "bias"  acceleration. 

In  the  diagram  of  Figure  II. 3,  we  have  a  complete 
representation  of  the  system  with  the  states  used  in 
modelling  our  missile. 

The  "STATE  EQUATIONS"  which  represent  our  system  are  the 
f ol lowing: 


x^=  -150.0  x^  -  2,646.66  x  -  705.75  x 


X  ,=  1 .  353  X  ^  +  X  ,  -  1 .  353  i)    +  1 .  334  H 

i  13  ZC.  S 


X  ,=  -6.572  x^  -  5  X,  +  6.572  i)    -6. 481  H 

3  13  'ZC  S 


X  =  -44.332  x^  -  59.109  x 

■  So 


x^=  x^  -  0.1482  x^  -  0.0395  x^  (II.  1) 

St  S  O 


x^=  -188.4  x^  +  188.4  x^ 
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x^=  -0.4608  X  ^-  2.231  x^  -  0.3406  x^  +  2.231  x   + 
-15.0949  x^  -  20.13  x^  -  0.143  x   +  0.4608  i) 
-  0.4544  H 


x.=  -50.0  x„  -  495.0  X   +  495.0  H 

S  8  9  5 


X  =  X 

9  10 


x^„=  -17.644  X   -  4.705  x 

10  S  fr 


This   system   o-f   equations  can  be  represented  in  matrix 
form,  where  we  have: 

X  =  A  X  +  B  u 


A  = 


-150  0     0 

1 . 35  0    1.0 

-6.5  O  -5.0 

0  0     0 

0  0     0 

0  0     0 

-0.46  -2.2  -3.4  2.2  -15.1 

0  0     0     0      O 

0  0     0     0      O 


O   2,646.66  -705.74   O 


O 
O 


O      0 
O      O 

O   -44.33    -59.11 
1.0   -0.15    -0.04 
O      O 


O    O 

O    O 

O    O 

O    0 

-188.4   188.4  O 

-20. 13   -O. 14  O 

O      O  -50. 

0      0    0 


0 
0 
O 
O 
O 
O 
0 


O 
O 
0 
O 
O 
0 
O 


17.64    -4.71 


O 


■495  O 
0  1 
0    O 
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F  = 


B 


0 
1.353 

1 

-6 

0 
.334 
.481 

0 

0 

0 

0 

0 

0 

0 

0 

0    0    0    0    0 
0    0    0    0    0 


0.461   -0.454 

O      495.0 

O       O 
u  =  -  F  X 

O     O     0.049  0.98615  1.404 
0     0       0        0       0 


The  vector  input  "u"  has  two  components 

D    —  the  commanded  acceleration:  and 
zc  ' 

H    -  the  desired  altitude. 

s 

The    input    d      is   considered   as   two   parts;   one 

zc  ^ 

corresponds    to   the   feedback   of   the   states   (from   the 

controller)  and  another  is  the  "bias"  acceleration  of  Ig. 
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TABLE  I 
DESCRIPTION  OF  THE  SYSTEM  STATES  AND 
AERODYNAMIC  COEFFICIENTS 

^^       filtered  acceleration 

X  2       acceleration  compensator 


X 


3 


X 


H 


acceleration  compensator 


q  -  pitch  angular  rate 


X5  on.   -    angle  o-f  attack 

X.  &  -  pitch  tail  incidence 

*  P 

x_  &   -  commanded  pitch  tail  incidence 

7  pC  ^ 

x^  X  -  control  1  er  compensator 

x^  H<t>  -altitude 

x^^  v(t)  -  velocity 


AERODYNAMIC  COEFFICIENTS 


C      slope  o-f  curve  o-f  pitching  moment  coef-ficient 
mot 

C    change  in  C   per  degree  pitch  control  incidence 
m  Sp  n\ 

C.,  slope  o-f  curve  o-f  normal  force  coefficient  C^.x  ot 
Noj  N 

C^,   change  in  C.,  per  degree  pitch  control  incidence 
Np      '       N        ^     "^ 
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The  model  can  be  separated  in  two  main  parts;  the 
CONTROLLER  and  the  AUTOPILOT/AIRFRAME. 

B.   THE  CONTROLLER 

The  controller  Mas  developed  based  in  the  procedure 
introduced  by  Dowdle  CRe-f.2D. 

The  missile  was  considered  as  a  point  mass  and, 
there-fore  the  altitude  as  -function  o-f  time  is  calculated  by 
double  integration  o-f  the  achieved  acceleration. 

The  controller  was  developed  without  take  in 
consideration  the  missile  itself,  i.e.,  the  trans-fer 
-function  G(s)  in  Figure  II. 2  was  assumed  as  "unity". 

According  to  the  CRe-f.2D,  this  procedure  is  valid  i-f  the 
bandwidth  o-f  G<S)  is  larger  than  the  crossover  -frequency  of 
the  controller. 

Due  to  the  constraints  of  the  scenario,  the  controller 
has  to  be  designed  in  such  way  that  no  overshoot  will  be 
tolerate.  Further  requirements  are  that  the  angle-of-attack 
and  the  acceleration,  have  to  remain  below  a  certain  level. 

The  controller  are  responsible  for  the  commanded 
acceleration  by  feedback  of  its  states. 

The  technique  used  to  solve  the  feedback  problem  was  the 
LINEAR  QUADRATIC  STATE  FEEDBACK  REGULATOR.  This  solution  has 
a  large  stability  margin  with  respect  to  phase  variations, 
as  well  infinite  gain  margin. 
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The  "cost  function"  to  be  minimized  is  the  following: 


:Tl 


J  =  [  <:  z  It)  Q  z  (t) 


*  uT 


t)  R  u(t)  >  dt     (II. 2) 


In   Figure   1 1. 4  we  have  the  schematic  representation  of 
the  CONTROLLER. 


H 


t) 

1 

V(t) 

1 

H(t) 

+ 

s 

s 

Figure  1 1. 4  -  The  Controller 


The  LEAD  compensator  indicated  in  the  controller  is  used 
in  order  to  improve  the  transient  response  characteristics. 
Its  transfer  function  is  : 


G  (s)   = 

c 


K*(s+w  ) 

z 

(s+w  ) 

p 


where 


K  =  10.0 


w  =  0.5  rad/sec 

z 


w  =  50.0  rad/sec 
P 


In   order  to  have  the  state  equations  of  the  CONTROLLER, 
we   draw   in   FIGURE   II. 5   its   signal  flow  graph  using  the 
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adequate   transformation   of  the  compensator  and  showing  the 
states  that  will  be  feedback. 


H 


D  (t)    1 


1/s    v(t)    1/s    H(t)    1 


'  1 


Figure  II. 5   Signal  Flow  Graph  of  the  Controller 


The  state  vector  is  defined  as: 


x(t)  =  C  x(t)   (H(t)-H  )   v(t)  1 


(II. 2) 


In   the   Figure   I I. 5   the   states  of  the  controller  ar« 
represented  as  indicated  in  Table  I. 

From  the  diagram,  we  have  the  following  STATE  EQUATIONSi 


X  =  A  X  +  B  u 


(II. 3) 


2  =  D  X 
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The  state  matrices  are; 


A   = 


P 

K(w    - 

z 

M      ) 

p 

0 

0 

0 

1 

0 

0 

0 

B   = 


O 
0 
1 


D   =    C  1  K  0  ] 
Since   we   consider   to   feedback  all  the  states,  let  in 
equation  1 1 . 1 

Q  =  I  (Identity) 


and 


R  =  0.04 


Using   the   OPTSYS  program  to  solve  that  OPTIMUM  CONTROL 
problem,  the  -following  -Feedback  gains  were  obtained; 

F   =  0.049 

F^  =  0.98615 

F   =  1.4044 

correspondents   to   the  feedback  of  the  states  indicated 
in  Figure  II. 5. 
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C.   THE  AUTOPILOT/AIRFRAME 

The  controller  Mas  developed  assuming  the  trans-fer 
•function  o-f  the  Autopilot/Air-frame  as  unity. This  assumption 
implies  that  the  achieved  acceleration  is  equal  to  the 
commanded  acceleration;  however,  they  are  di-f-ferent  due  to 
the  aerodynamic  of  the  missile. In  order  to  take  that  in 
consideration  we  have  to  introduce  the  autopi lot/airf rame 
model . 

Normally,  the  sea-skimming  is  a  SKID-TO-TURN  missile 
(turns  instantaneously)  ,  but  due  to  problems  o-f 
classification  we  have  used  for  our  Sea — Skimming  model  the 
Autopilot  given  by  Arrows  CRef.lD  for  a  BANK-TO — TURN 
missile  (has  to  bank  in  order  to  turn)  that  was  designed  to 
fly  at  30,000.0  ft  instead  of  sea  level. 

Take  in  consideration  that  the  main  objective  of  this 
work  is  to  establish  a  procedure  for  sensitivity  and 
robustness  analysis  and  that  we  are  studying  only  the 
vertical  movement  of  the  missile  uncoupled  of  any  other 
movement, we  consider  that  the  approximation  is  useful. 

In  Figure  II. 6  we  have  the  block  diagram  of  the  missile, 
representing  the  plant  without  the  controller. 

The  state  equations  derived  from  this  diagram  are  on  the 
Appendix  A. 
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Figure  II. 6   Block  Diagram  of  the  Airf rame/Autopi 1  at 
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D.   THE  TIME  RESPONSE 

Using  the  so-Ftware  CONTROLS,  the  system  represented  by 
equations  (II. 1)  Mas  simulated  and  the  time  response  for  the 
altitude,  velocity,  acceleration  and  angle  of  attack  Mere 
obtained. 

The  altitude  is  in  Figure  II. 7  and  Me  can  see  that  the 
missile  arrives  at  the  desired  altitude  in  about  10  seconds. 

No  overshoot  is  shOMn  and  a  small  steady-state  error  is 
presented  characterizing  our  Type  0  system. The  steady-state 
error  Mas  less  than  5/C. 

The  angle-of -attack  , reproduced  in  Figure  1 1. 8,  varies 
betMeen  3.7  deg  and  -1.2  deg  satifying  the  condition 
required  by  the  airframe,  as  indicated  in  CRefs.  1  and  2D. 

The  acceleration  is  in  Figure  II. 9  and  the  velocity  is 
in  Figure  11.10  representing  a  smooth  behavior  of  the 
system. 

In  the  next  chapter  Mill  be  presented  the  Sensitivity 
Analysis  using  the  developed  model. 
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III.  SENSITIVITY  ANALYSIS 

A.    INTRODUCTION 

The  objective  o-f  the  SENSITIVITY  ANALYSIS  is  to 
analyse  how  the  behavior  o-f  the  system  is  af-fected  by 
variations  o-f  the  aerodynamic  parameters.  Those  variations 
are  due  to  uncertainty  of  the  design  and/or  problems  ocurred 
during  the  flight. 

From   the   developed  model  we  can  see  that  the  important 

parameters   to   our   system   are   C   ,  C^.  ,  C    and  C^,  ,   as 
^  '  moj'   Not'   m&       N&* 

described  on  TABLE  I. 

The  variations  could  be  interpreted  as  perturbations  in 
our  STATE  EQUATION  SYSTEM  (A-matrix)  and  the  analysis, 
basically,  take  into  consideration  the  change  in  the  poles 
positions  due  to  change  in  the  parameters. 

As   we   have   feedback,   we   replace  the  A  matrix  by  the 

augmented  A-matrix (  A   =  A  -  B  F  ). 
^  a 

The  sensitivity  of  a  matrix  where  the  elements  are 
time-invariant , can  be  verified  in  different  ways.  Among 
those  procedures,  some  are  enumerated  below. 

First,  using  the  method  introduced  by  Golub  and  Van  Loan 
CRef.33  we  can  analyse  the  matrix  to  verify  how  it  is  "ill- 
conditioned";  where  an  "ill-conditioned"  matrix  means  a  very 
sensitive  system. 
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In  order  to  accomplish  the  analysis,  we  have  to 
calculate  the  "condition  number".  It  is  defined  as: 

»<(A  )  =  II  All  II  A~^|  I 
a  a         a 

where 

II  A^  I  I  is  the  norm  o-f  augmented  A  and 

II  A   II  is  the  norm  o-f  inverse  augmented  A 

The  condition  mumber  depends  upon  the  considered  norm 
but  it  can  be  proven  that  the  condition  of  the  matrix  will 
be  the  same  for  all  norms. 

When  we  use  the  2-norm  one  has 

K  (A  )  =  "    ^^a^  (III.  1) 

s     (A  ) 

-    a 

where  a  (A  )  and  £  (A  )  are  the  maximum  and  minimum 
singular  value  of  the  A-augmented  matrix,  respectively. 

Calculating  the  condition  number  for  the  different 
values  of  the  parameters  it  is  possible  to  establish  how  the 
sensitivity  of  the  system  changes  when  there  are  changes  in 
some  elements  of  the  state  matrix. 

The  problem  of  this  kind  of  analysis  is  that  we  can  not 
verify  the  sensitivity  of  the  system  with  respect  to  a 
specified  element  of  the  state  matrix  but  only  the 
sensitivity  of  the  overall  system. 

In  Frank  CRef.4],  we  have  another  method  where  the 
analysis   is  done  based  upon  the  STATE  SENSITIVITY  EQUATIONS 
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that  could  be  characterized  in  time  domain,  frequency  domain 
orin  terms  of  a  performance  index. 

The  basic  idea  is  to  define  "sensitivity  functions"  and 
analyse  the  behavior  of  the  system  by  those  functions. 

The  analysis  starts  with  the  system  defined  by: 

X  =  f  (x,t,u,oc)  (111.2) 

where: 

-  x_   represents  the  state  vector; 

-  t  is  time  ; 

-  u  is  the  input  vector  ;  and 

-  ai   represents  a  vactor  of  the  different  parameters. 

If  o^  is  changed  by  Aoc  (some  perturbation)  and  using 
Taylor's  expansion, we  have: 


,«)  =2_ 


Ax  (t,o<)  =)  Ax   I   ioc.  (III. 3) 

The   subscript   q<„  indicates  that  the  partial  derivative 

is  taken  at  the  nominal  parameter  values. 

The  partial  derivative  is  called  the  trajectory  vector  ( 

X   )   with   respect   to   the   j —  parameter  and  has  the  same 

dimension  as  the  state  vector. 

The   components   of   that   vector  are      defined   as   the 
TRAJECTORY  SENSITIVITY  FUNCTION. 


X.  .(t,oc  )  =   t^^^^°^^   I  (III. 4) 
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The   equation  III. 4  represents  the  partial  derivative  o-f 

.  .    . th   .  .  th 

the  1 —  state  with  respect  to  the  j —  parameter.  Assuming  we 

have   n  states  and  m  parameters,  it  will  be  possible  to  -form 

a  n  X  m  matrix  called  TRAJECTORY  SENSITIVITY  MATRIX  (III.5). 

The  matrix  is  also  called  the  JACOBIAN  MATRIX. 


^xn 
doc 


ax  1 
m 


axn 

doc 
m 


(III. 5) 


The      change       in      the      behavior      o-f       the      system      will    be 
determined    by: 


Ax    =    X     (t,o«„)     Aoc   =^     X.    Aoc . 
-        -  -»         -      Z_.      J         J 


(III. 6) 


J  =  i 


The  analysis  is  based  in  how  the  sensitivity  function 
changes  when  we  change  the  parameters.  The  result  could  be 
visualized  by  plotes  o-f  X  as  -function  o-f  time,  as  states  of 
the  expanded  system. 

Finally,  we  consider  the  EIGENVALUE  SENSITIVITY 
presented  by  Porter  &  Crossley  CRef.53. 

This  latter  procedure  was  choosen  to  be  used  in  the 
analysis  of  this  thesis,  where  the  sensitivity  of  the  sistem 
is   related   to   the   variation   of   the   eigenvalues  of  the 
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augmented   matrix   (A    =   A   -   B   F)   with   respect  to  the 

variation  o-f  the  parameters. 

In  order  to  -further  develop  our  analysis,  the  concept  o-f 
MODAL  ANALYSIS  will  be  introduced. 

B.    THE  MODAL  ANALYSIS 

The  concept  is  based  in  generating  the  input  vector  o-f  a 
system  by  linear  -feedback  of  the  state  vector  in  such  way 
that  the  prescribed  eigenvalues  are  associated  with  the 
dynamic  model  o-f  the  resulting  close— loop  system. 

As  an  example,  let  us  assume  a  scalar  system: 

X  =  a  X (t)  +  b  u(t) ;  (III. 7) 

i-f  u(t)  =  O  ====>  X  =  X   exp  (at) 

where  exp  (at)  de-fines  the  "mode"  o-f  the  system. 

If  a  ^  O  the  system  is  stable  and  if  a  ^  O  the  system  is 
unstable. 

Assuming  g  as  the  feedback  gain,  we  have: 

u(t)  =  g  X (t)  (III. 8) 

Inserting  the  value  of  u  from  equation  1 1 1. 8  in  equation 
III. 7: 

X  =  (a  +  b  g)  X  (t)  (III. 9) 

and,  solving  for  x(t): 

x(t)  =  X (O)  expC(a  +  bg)tD 
Now  the  mode  of  the  system  is  defined  by  the  exp(a+bg). 
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When  Me  have  a  state  vector,  the  system  is  represented 
by 

X  =  A  X  +  B  u  (III. 10) 

Considering 

u  =  -  F  X 
and    assuming    all    the    states   as   observable   and 
controllable,   the   substituting   value   o-f   u   in   equation 
(III. 10) ,  Me  have 

x=(A-BF)x    or    x=A   x        (III. 11) 

a 

The  precise  nature  of  the  -free  motion  of  the  continuous 
time  system  folloMing  any  disturbance  can  be  described  in 
terms  of  eigenvalues  and  eigenvectors  of  the  augmented  plant 

matrix  (A  ) . 

a 

Assuming  A   has  n  distinct  eigenvalues  (x  ,x  ,...,x  ); 

then   it   also   has   n   corresponding   linearly   independent 

eigenvectors  u  ,u  ,...,u   such  that 
12       n 

A   u.  =  X.  u. ;  i=l,2,...,n  (III. 12) 

a   1     1   1 '     '  '    ' 

Using  the  calculated  eigenvalues  and  eigenvectors  Me  can 

Mrite  the  modal  matrix  of  A   as: 

a 

U  =  Cu   u  .  .  .    u    1 
12       n 

Mhere  each  column  of  U  is  the  eigenvectors  of  A  . 

^  a 

Based  in  the  modal  matrix  U,  the  equation  (III. 12)  could 

be  Mritten  as 
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A   U  =  A  U  (III.  13) 

a 

with  A  as  a  diagonal  matrix  o-f  the  eigenvalues. 
The  matrix  A  is  obtained  from  equation  (III. 13): 

U~^  A   U   =   A  (III. 14) 

a 

In    addition    to   the   eigenproperties   o-f   A  ,   the 
corresponding   properties  o-F  its  transposed  matrix  (A')  play 

CI 

an  important  role  in  the  modal  analysis. 

The   matrix   A'   has   the   same  eigenvalues  of  A  ,  but  a 
a  ^  a 

di-f-ferent  sets  of  eigenvectors.  It  can  be  represented  as 

A'  V.  =  X.  v.:   j  =  1,2,...  ,n 

where  v.  represents  the  eigenvectors  of  A'. 
J  ^  a 

In  CRef.5D  it  is  shown  that 

V'.  u.  =  u.'  V  .  =  &.  .  (III.  16) 

J   X      1   J      ij 

where  &.  .  is  the  Kronecker  delta:  i.e. 

&.  .  =  1       if   i  =  j 

X  J 

and         &.  .  =  0       if   i  9«  j 
X  J 

Thus  in  matrix  form 

V'U  =  I  (III. 17) 

where  the  collumns  of  V  are  the  vectors  v.. 

1 

Using  the  modal  matrix  of  A   and  introducing  a  new  state 

a  ^ 

vector  y^(t)  ,  we  have: 

X (t)  =  U  y (t)  (III. 18) 
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The  derivative  of  both  sides,  gives  us: 
x<t)  =  U  y(t) 


with 


X  (t)  =  A  X  <t) 

a 


(111. 19) 

(111. 20) 


Combining   equations   ( 1 1 1 . 18) , ( 1 1 1 . 19)  and  (III. 20),  Me 
have: 


—  1 


y  =  U    A   U  y(t) 


and  using  the  equation  (III. 14), 


y(t)  =  A  y(t)  (III. 21) 

Because   the   matrix   A   is  diagonal   the   equations   are 

uncoupled   and   the   solutions  of  the  system  represented  by 

(III. 16)  are  given  by 

y. (t)  =  y. (O)  exp(x.t)  (III. 22) 

1        1  1 

Aplying  equation  III. 18, 


X.  (t)  =  Cu,  u, 

1  12 


.  u  : 

n 


n 


(II 1.23) 


which  implies  that 

x(t)  =  u,y, (0)  exp(X,t)  +  u„y„(0)  exp(X  t)  + 


I'l 


2'2 


Taking  the  value  of  y(t)  in  equation  (III. 18), 


y(t)  =  U  ^x (t) 


Using  equation  (III. 17), 


y.  (O)  =  V.'  X  (o) 


(III.24) 


(III. 25) 
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Combining  equations  (I  II. 23)  and  (II  I. 25): 


x(t)  =  \  X (O)  u.v:  exp(X.t)  (III. 26) 

/        1  1    '^   1 


Equation   1 1 1 . 26   shoNS   that   the  system  represented  by 

(I  II. 8)   has   its   modes  described  by  the  eigenvectors  of  A 

a 

and   its   transpose   and   eigenvalues  of  A   matrix.  Also,  in 

order   to   have  the  system  stable,  x(t) >  O  as  t >  a>  , 

which  happens  if  and  only  if 

Re  X.  <  O  ,  for  all  the  values  of  i 

In  the  Modal  Control  CRef.5]  all  this  procedure  is  shown 
for  the  case  when  some  of  the  eigenvalues  are  equal  which 
does  not  apply  to  our  system. 

As  the  role  of  the  eigenvalues  is  so  important,  we  can 
verify  how  the  sensitivity  of  our  control  1  system  by 
analysing  the  change  in  the  eigenvalues  of  the  matrix  A  or 
the  augmented  A  for  the  close  loop  system  when  any  element 
is  changed. 

In  order  to  avoid  calculating  the  new  eigenvalues  for 
each  different  A,  the  equation  (III. 26)  can  be  used  to  find 
approximated  values  that  can  be  analysed  analytically. 
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C.    EIGENVALUE  SENSITIVITY 

In  our  system  we  consider  the  change  in  the  elements  of 
the  matrix  due  to  variations  in  the  aerodynamic 
coe-f -f  icients. 

Basically,  we  concerned  with  the  eigenvalues  that  are 
changed  in  the  direction  of  the  right  hand  plane, making  the 
system  unstable  or  affecting  the  time  response. 

The  result  of  the  analysis  could  be  plots  or  tables 
relating  the  eigenvalues  with  the  variation  of  the 
parameters. 

For  each  change  in  the  parameters,  we  have  a  new  matrix. 

Approximations  to  the  eigenvalues  associated  with  the 
different  plant  matrices  in  the  vicinity  of  the  original  one 
could  be  used  without  having  to  compute  the  eigenvalues  of 
each  new  matrix. 

Me  describe  the  approach  presented  in  the  reference 
CRef.5]  applying  the  "first  order  eigenvalue  sensitivity"  to 
the  change  of  one  of  the  parameters  and  doing  the  complete 
analysis  with  the  calculated  eigenvalues. 

1.   First-order  Eigenvalue  Sensitivity 
From  equation  III.  11,  we  have: 

x(t)  =  A^  x(t) 


Let 's  assume 


A 


=  [  %^] 


where  a,  ,  is  any  element  of  A^. 
kl  a 
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Di-f f erentiating   equation   (III.  12)  with  respect  to 


a. ,  indicates  that 
kl 


dfii      u.     +    A    aui       =    axi    u.     -»■    X       aui  (1 11.27 

da,  ,  da,  .  da,  .  da,  . 

kl  kl  kl  kl 

multiplying  (III. 27)  by  v.'  then  gives 

V .'  dA  u .  +  V  f  A  dui   =  V .'  dxi  u .  +  X .  v .'  dui 
1  T—   1     1 1 1     11 

^^kl  ^^kl        ^^kl  ^^kl 

As    V.'    A    is    equal    to    x .     v.'     ,    we    have 
1  ^  I       1     ' 

V.'     dA    u.     =    V.'     dxi    u. 
1    -r—      I  1    T^ 1 

^^kl  ^^kl 

Which  reduces  to  the  set  o-f  equations 

dxi  =  V. (k)  u.(l)   ;  i=k=l=l,2,. . . ,n   (III. 28) 
T —     1       1  1     y  1 

^^kl 

Due  to 

dA   =   d.,  d., 
^^kl 


In  equation  (III. 28),  v.(k)  and  u.(l)  represent  the 

k —  element  of  v.  and  the  1 —  element  of  u. ,  respectively. 

Those    coefficients    may   be   considered   as   the 

elements  of  a  set  of  n  eigenvalue  sensitivity  matrices. 

The   new   eigenvalues   (x.)   of  the  altered  matrix, 

1 

when  the  element  a.,  is  changed  will  be  calculated  as 

kl 

^  d  X 

\  =  S  "  — i  ^  ^ki 

^^kl 
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When  we  have  more  than  one  parameter  altered,  the 
linearity  could  be  used  and  the  final  eigenvalue  Mill 
obtained  -from  the  summation  of  the  sensitivity  coefficients 
multiplied  by  the  respective  change  in  the  coefficients. 

As  an  example  of  this  procedure,  let's  consider  our 

system,  changing  the  value  of  C    such  that 

msp 

AC  =0.02 

m&p 

Mhen      C  is      changed      as    specified,    the    fol lowing 

m&p 

elements      of       the      augmented-A    matrix    of    the    state    equations 

presented    in    the    Chapter    II    are   changed: 

A    (4,6) :  Aa„^    =    14.778 

A     (7,6) :  Aa,,    =      5.035 

a  76 

The  eigenvalues  and  eigenvectors  of  the  augmented  A 

matrix  and  of  its  tranpose  were  calculated  using  the  program 

Mat lab. 

a)  Eigenvalue  X   : 


X   =  -159.75  +  J  18.85 


A  X   =Cv  (4)  u  (6):  A  a„^   +  Cv.(7)  u  (6) 3  A  a 


V  (4)  =  -0.0148  -  j  0.0039 


v^(7)  =  1 


u  (6)  =  0.014  -  j  0.027 


X      =  -159.815  +  j  18.85 
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b)  Eigenvalue  X_  : 


=  -159.815  -  j  18.85  (complex  conjugate  o-f  X  ) 


c)  Eigenvalue  x   : 


Xj   =  -50.0 


AX,  =  Cv,(4)  u,(6)]  Aa„^  +  Cv,(7>  u,(6)D  Aa,^ 

3  3         3            Ho         3         3            To 

V  (4)  =  0.0003 

Vj(7>  =  0.0039 

Uj(6)  =  -0.001 


X^   =  -49.99998 


d)  Eigenvalue  X 


X^  =  -8.61  +  j  7.91 

AX^  =  Cv^(4)  u^<6)D  Aa^^  +  Cv^  (7)  u^  (6>  ]  Aa-,^ 

v^(4)  =  -0.1667  -  j  0.1322 

v^(7)  =  1 

u^(6)  =  -0.0334  -  j  0.132 

X„     =  -8.95  +  j  7.6355 


e)  Eigenvalue  X   : 


X,     =  -8.95  -  j  7.6355 


f)  Eigenvalue  X 
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X^  =  -1.82  +  j  2.02 

AX  =  Cv^<4)  u^<6)D  Aa    +  Cv  (7)  u  (6) D  Aa 

v^(4)  =  0.1718  -  j  0.0161 

v^(7>  =  -0.5737  +  j  0.1689 

u^(6)  =  0.01  +  j  0.0363 


X^     =  -1.846  +  j  2.05 


g)  Eigenvalue  x 


X^     =  -1.846  -  j  2.05 


h)  Eigenvalue  X, 


X.  =    -2.71 

8 


AX„         =    Cv„(4)    u^(6):    Aa         +    Cv    (7)    u    (6)3    Aa 

8  88  Hfr88  76 

V  <4)      =     O. 1741 

8 

V  (7)     =    -0.6597 

8 

u.(6)     =    -0.0396 

8 


X.  =    -2.681 

8 


i)  Eigenvalue  x.  : 

^  9 


X.     =  -0.5 

9 

AX     =  Cv  (4)  u.(6):  Aa„^  +  Cv„(7)  u„(6)]  Aa,^ 

9  9  9  Hfr         9  9  76 

V  (4)  =  0.1718  +  j  0.016 

9 
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V  (7)  =  -0.5737  -  j  0.1689 

9 


u  (6)  =  0.01  -  j  0.0363 

9 


X„     =  -0.4741 

9 


j)  Eigenvalue  x  ^  : 


So    =-0.14 

^10    =  ^^iO<^^   ^10^^^^  ^H.  -^   '=^iO<^>   ^10<^>^   ^7. 

V  (4)=  O. 1741 

1  0 

V  (7)=  -0.6597 

1  0 

u   (6)=  -0.0396 

1  0 


X, „    =  -O. 1104 
1  0 

For  the  analysed  parameter,  our  system  is 
unsensitive,  taking  in  consideration  the  small  di-f -f erences 
between  the  new  and  original  eigenvalues. 

Repeating  this  procedure  for  all  different  parameter 
values  Mill  give  us  the  approximation  of  the  new  eigenvalues 
and  that  can  be  helpfull  in  a  qualitative  analysis  of  what 
parameter  has  more  effect  over  our  model . 

As  we  want  to  verify  the  effect  in  the  time 
response  and  due  to  the  Software  available  the  analysis  will 
be  done  completely  with  the  calculated  exact  values. 
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2.   Exact  Eigenvalues 

Using  the  Controls  Program,  we  will  veri-fy  the  variation 
o-f  the  eigenvalues  of  the  augmented  A  matrix  (  A  -  FB)  with 
respect  to  the  variation  o-f  the  parameters  and  by  simulation 
we  analyse  the  e-f-fects  on  the  time  response. 

At  the  end  o-f  the  analysis  we  compare  the  results  oi  the 
example  used  for  the  first  order  approximation  with  the 
exact  values  that  were  calculated. 

As  reference,  the  eigenvalues  of  the  original  system 
are: 


X   =  -159.746  +  j  18.979 


X   =  -159.746  -  j  18.979 

2 


X^  =  -49.9992 


X^  =  -8.60697  +  j  7.90968 


X   =  -8.60697  -  j  7.90968 


X   =  -1.81854  +  j  2.0238 


X   =  -1.81854  -  j  2.0238 


X   =  -2.70266 

8 


X   =  -0.503819 

9 


X   =  -O. 143073 

1  0 
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Each   aerodynamic   coe-f  f  icient ,   that   takes  part  of  the 
considered   model  was  changed  up  to  +25%  and  the  eigenvalues 

o-f   the   perturbed  system  matrix  < A— augmented)  were  analysed 
as  well  as  the  TIME  RESPONSE  of  the  system  . 

a.  C 

mo( 

On   TABLE   II  we  have  the  eigenvalues  related  to  the 

considered    values   of  C   and  as  we  can  see,  with  exception 

moc  '  ^ 

of   the  7 —  eigenvalue,  the  changes  are  small  and  the  system 
is  pratically  insensitive  to  that  variation. 

The  time  responses  are  shown  in  Figures  III.l,  111.2,111.3 
and  II 1.4  . 

-  Angle-of-attack   (Figure   III.l)   -  changes  in  the  same 

direction   of  C    with  the  minimum  value  changing  about 

moc  -s   -s 

0.5   degrees   in   the   negative   direction.  The  maximum 
value  stays  almost  constant  ; 

-  Acceleration   (Figure  III. 2)  -  the  change  was  less  than 

2 

5  ft/sec   either  the  maximum  or  the  minimum  values; 

-  Velocity  (Figure  I I I. 3)  -  almost  the  same  ;  and 

-  Altitude  (Figure  III. 4)  -  almost  the  same. 

b.  C 

m&p 

The   change   in   the   eigenvalues  are    on  TABLE  III  ;  the 
time  response  is  on  Figure  111.5,111.6,111.7  and  III. 8. 

The   variation   of   the   eigenvalues  are  slightly  higher 

with  changes  in  C    but  the  system  stays  stable  and  the  time 

^        msp  ^  ^ 

response  has  a  small  variation. 
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-  Angle-of-Attack  (Figure  III. 5)  -  The  minimum  values 
increases  negativelly  about  0.5  degrees  but  in  the 
contrary  direction  o-f  the  coe-f  f  icient ,  i.e.,  increasing 
the  value  o-f  C^    negativelly,  the  minimum  Mas  smaller; 

-  Acceleration  (Figure  III. 6)  -  small  changes; 

-  Velocity  (Figure  III. 7)  -  small  changes; 
Altitude  (Figure  III. 8)  -  pratically  the  same. 

The  change  in  the  eigenvalues  are  on  TABLE  IV  and  the 
time  response  is  on  Figures  111.9,111.10,111.11  and  III. 12. 

The  changes  in  the  eigenvalues  are  small  but  the  e-F-fect 
on  the  time  response  is  more  sensible,  mainly  on  the  angle— 
o-f —attack. 

-  Angl  e-o-f-Attack  (Figure  III. 9)  -  the  maximum  and  the 
minimum  increase  1.0  degree  in  the  positive  and 
negative  directions,  respecti vel ly; 

-  Acceleration  (Figure  III.  10)  -  maximum  changes  o-f  10/C 
in  the  negative  direction  (minimum)  and  20'/.  in  the 
positive  direction  (maximum) . These  values  are  smaller 
when  we  increase  the  value  o-f  C.,   ; 

-  Velocity  (Figure  III.  11)  -  small  changes; 

-  Altitude  (Figure  III. 12)  -  does  not  change. 

N&p 

The  change  in  the  eigenvalues  are  on  Table  V  and  the 
time  response  Are  on  the  Figures  111.13,111.14,111.15  and 
III. 16 

-  Angl  e-o-f-Attack  (Figure  III.  13)  -  Changes  almost  1.0 
degree  with  the  maximum  and  minimum  increasing  when  the 
value  o-f  the  parameter  increases; 
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-  Acceleration   (Figure   III. 14)   -  the  maximum  increases 
Mith  the  value  o-f  the  parameter; 

-  Velocity   (Figure  III.  15)  -  pratically  the  same  -for  all 
the  considered  values  o-f  the  parameter; 

-  Altitude  (Figure  III. 16)  -  constant. 
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TABLE  II 

C    -  SENSITIVITY 

C           -0 

.04 

-0.  05 

-159.741 

+ 

J 

18.9843 

-159.743 

+  j  18.9816 

-159.741 

- 

J 

18.9843 

-159.743 

-  j  18.9816 

-49.9992 

-49.9992 

-8.55822 

+ 

J 

6.37130 

-8.56386 

+  j  7. 18736 

-8.55822 

- 

J 

6.37130 

-8.56386 

-  j  7. 18736 

-2.31740 

+ 

j 

2.95574 

-2. 14117 

+  j  2.43275 

-2.31740 

- 

j 

2.95574 

-2. 14117 

-  j  2.43275 

-1.79944 

-2. 14200 

-0.51650 

-0.510282 

-0. 14304 

-0. 14306 

C 
m 

-0 
at 

.07 

-0.08 

-159.748 

+ 

J 

18.9762 

-159.751 

+  j  18.9735 

-159.748 

- 

j 

18.9762 

-159.751 

-  j  18.9735 

-49.9992 

-49.9992 

-8.66197 

+ 

J 

8.55389 

-8.72178 

+  j  9. 14021 

-8.66197 

+ 

j 

8.55389 

-8.72178 

+  j  9. 14021 

-1.47981 

+ 

j 

1.81088 

-1.22077 

+  j  1.70081 

-1.47981 

- 

j 

1.81088 

-1.22077 

-  j  1.70081 

-3.27131 

-3.67125 

-0.49779 

-0.49148 

1 

-0. 14308 

-0. 14308 

51 


TABLE 

III 

C    -  SENSITIVITY 
m&p 

c 

m&p 

-0.06 

-0.07 

-162.604 

+  j  16.7131 

-161.200   + 

j  17.9141 

-162.604 

-  j  16.7131 

-161.200   - 

j  17.9141 

-49.9992 

-49.9992 

-5.74603 

-»-  j  8.86120 

-7.13033   + 

j  8.50992 

-5.74603 

-  j  8.86120 

-7. 13033   - 

j  8.50992 

-1.24694 

+  j  1.66681 

-1.54230   + 

j  1.80335 

-1.24694 

-  j  1.66681 

-1.54230   - 

j  1.80335 

-3.86370 

-3.30548 

-0.49103 

-0.49828 

-0. 14311 

-0. 14309 

C  ^ 

-0.09 

-0.  10 

-158.243 

+    j     19.9192 

-156.678   -»- 

j  20.7492 

-158.243 

-  j  19.9192 

-156.678   - 

j  20.7492 

-49.9991 

-49.9990 

-10. 1621 

+  j  6.98960 

-11.7920   + 

j  5.53519 

-10. 1621 

-  j  6.98960 

-11.7920   - 

j  5.53519 

-1.47981 

+  j  1.81088 

-2.03752   + 

j  2.55411 

-1.47981 

-  j  1.81088 

-2.03752   - 

j  2.55411 

-2.27254 

-2.02276 

-0.50786 

-0.51120 

-0. 14307 

-0. 14306 
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TABLE 

IV 

C^^   -    SENSITIVITY 

•^Noc 

0. 

13 

0.  14 

-159.750 

+ 

J 

18.9895 

-159.748   +  j 

18.9842 

-159.750 

- 

J 

18.9895 

-159.748   -  j 

18.9842 

-49.9992 

-49.9992 

-8.76667 

+ 

J 

8.27012 

-8.68784   -•-  j 

8.09494 

-8.76667 

- 

J 

8.27012 

-8.68784   -  j 

8.09494 

-1. 18629 

-«■ 

j 

1.74321 

-1.47163   +  j 

1 . 84903 

-1.18629 

- 

j 

1.74321 

-1.47163   -  j 

1 . 84903 

-3.63456 

-3.22775 

-0.48878 

-0.49663 

-0. 14314 

-0. 14310 

'^Not 

0. 

16 

0.  17 

-159.744 

+ 

j 

18.9736 

-159.741  +    j 

18.9683 

-159.744 

- 

J 

18.9736 

-159.741   -  j 

18.9683 

-49.9992 

-49.9992 

-8.52441 

+ 

J 

7.71274 

-8.44077   +  j 

7.50213 

-8.52441 

- 

j 

7.71274 

-8.44077   -  j 

7.50213 

-2. 15317 

+ 

j 

2.32249 

-2.40105   +  j 

2.67918 

-2. 15317 

- 

j 

2.32249 

-2.40105   -  j 

2.67918 

-2.20610 

-1.88565 

-0.51041 

-0.51651 

-0. 14305 

-0. 14303 
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TABLE  V 

- 

C^,   -  SENSITIVITY 
N&p 

^N&p 

0.02 

0. 

03 

-159.004 

+  j  12. 1138 

-159.748 

+ 

j 

18.9842 

-159.004 

-  j  12. 1138 

-159.748 

- 

j 

18.9842 

-49.9996 

-49.9992 

-9.37545 

+  j  7. 11519 

-8.68784 

+ 

j 

8.09494 

-9.37545 

-  j  7. 11519 

-8.68784 

- 

J 

8.09494 

-2.24432 

+  j  2.56109 

-1.47163 

+ 

j 

1 . 84903 

-2.24432 

-  j  2.56109 

-1.47163 

- 

j 

1 . 84903 

-1.98723 

-3.22775 

-0.51356 

-0.49663 

-0. 14304 

-0. 14310 

^N&p 

0.05 

0. 

06 

-160. 106 

+  j  21.5632 

-160.459 

+ 

J 

23.8437 

-160. 106 

-  j  21.5632 

-160.459 

- 

j 

23.8437 

-49.9990 

-49.9988 

-8.25204 

+  j  8.25703 

-7.91449 

+ 

j 

8.57438 

-8.25204 

-  j  8.25703 

-7.91449 

- 

j 

8.57438 

-1.57370 

+  j  1.85345 

-1.35009 

+ 

J 

1.75196 

-1.57370 

-  j  1.85345 

-1.35009 

- 

j 

1.75196 

-3. 18783 

-3.60970 

-0.49848 

-0.49275 

-0. 14310 

-0. 14312 
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D.   CONCLUSION 

As   a   conclusion   for  the  sensitivity  analysis,  we  have 
the  -Following: 

With   respect   to   the   eigenvalues,   the   pitch  control 

system   is   more   sensitive  to  changes  in  C     and  C     than 

m  &p       N  &p 

changes  in  C    and  C^.    :  and 
man  Noc   ' 

Considering   the   time  response,  we  can  verify  that   the 

angle-o-F-attack   is   sensitive   to   all   the  parameters  with 

higher    variation    with   C^,  :   the   acceleration   is   more 

sensitive   to   C     and   C^,    ;  the  velocity  has  the  maximum 

mot         Noc  ' 

value   changed   less   than   107.   with   respect   to   all   the 

parameters;    and    the    altitude   is   insensible   for   the 

considered  variations. 

It   is  not  necessary  to  emphasize  the  importance  of  this 

analysis   for  the  design  of  the  control  system,  but  thinking 

only   in   the   eigenvalues   and  comparing  the  results  of  the 

first   order   approximation  presented  for  C  ,   with  those  of 

'^^  '^  m  &p 

Table  III  for  C    changed  to  -0.06  (AC    =  0.02)  we  can  have 

m  &p  m  &p 

an   idea   of   the   importance   of   the  parameter,  using  that 

approximation. 
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Figure  I I I . 1 


Angle-of  Attack  vs  Time 
■for  different  values  of  C 
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Figure  III. 2 


Acceleration  vs  Time 
for  different  values  of 
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Figure  I I I . 3 


Velocity  vs  Time 

■for  different  values  of  C 
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Figure  III. 4 


Altitude  vs  Time 

for  di-F-ferent  values  of  C 
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Figure    III. 5 


Angle-Q-f-Attack    vs    Time 
for    different    values    of    C 
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Figure  III. 6   Acceleration  vs  Time 

for  different  values  of  C 
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Figure    III. 7      Velocity    vs    Time 

for    di-f-ferent    values    o-f    C 
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Figure    III. 8 


Altitude    vs    Time 

for    di-fferent    values   of    C 
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Figure  III. 9   Angle-o-f-Attack  vs  Time 

■for  different  values  of  C 
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Figure  III. 10   Acceleration  vs  Time 

for  different  values  of  C 
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Figure  III. 11   Velocity  vs  Time 

for  different  values  of  C 
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Figure  III.  12 


Altitude  vs  Time 

for  different  values  of  C 
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Figure  III.  13   Angl e-o-f-Attack  vs  Time 

•far  different  values  of  C 
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Figure  III. 14   Acceleration  vs  Time 

•for  different  values  of  C 
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Figure  III. 15   Velocity  vs  Time 

•for  different  values  of  C 
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Figure    III. 16       Altitude    vs    Time 

-for    different    values    of    C 
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IV.  ROBUSTNESS  ANALYSIS 

A.   INTRODUCTION 

In  the  previous  chapter  we  take  into  consideration  the 
change  in  the  eigenvalues  and  in  the  time  response  o-f  the 
system  due  to  variations  on  the  aerodynamic  parameters  that 
take  part  into  the  model;  under  the  robustness  analysis  we 
veri-fy  what  we  can  expect  o-f  the  system  when  it  is  a-f-fected 
by  any  kind  o-f  perturbation. 

The  main  idea  is  to  use  the  robustness  analysis  to 
design  a  new  set  o-f  -feedback  gains  in  order  to  have  the  same 
behavior  in  di-fferent  situations  and/or  environments  which 
act  as  perturbations  to  our  missile. 

In  classical  -frequency  domain,  -for  single-input  single- 
output  systems (SISO) ,  a  robust  design  can  be  achieved  using 
Bode,  Nyquist  or  Nichols  plots.  With  these  techniques  it  is 
possible  to  de-fine  gains  that  give  us  gain  and  phase  margins 
for  a  "robust"  system.  For  multi-input  multi-output  systems 
(MIMO) ,  the  classical  techniques  are  no  longer  valid. 

Taking  this  in  consideration,  we  have  analysed  the 
robustness  o-f  the  original  design  by  means  o-f  the  minimum 
singular  value  o-f  the  return  di-fference  matrix  -for  di-f-ferent 
•frequencies  and  we  will  use  the  same  technique  to  improve 
the  control  system. 
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In  order  to  show  why  we  are  applying  this  technique,  we 
start  with  the  Nyquist  Criterion  -for  a  SISO  system  and 
extend  the  analysis  to  a  MIMQ  system,  using  the  minimum 
singular  value  theory,  as  explained  by  Lehtomaki ,  Sandel 1 
and  Athans  in  CRe-f.6]. 

B.   THE  NYQUIST  CRITERION 

Let's  take  in  consideration  the  SISO  system  represented 
in  Figure  IV. 1 . 


U 


■^Or 


Figure  IV. 1   Single-Input  Single-Output  System 

where  G(s)  includes  the  controller  and  plant  dynamics 
with  unity  -feedback  . 

The  Nyquist  criterion  states  that  if  the  open  loop 
trans-fer  -function  G(s)  does  not  have  any  pole  in  the  right 
hal-f  s  plane,  then  the  locus  o-f  G(s)  will  not  encircle  the 
point  (-1,0)  in  the  Nyquist  plot,  where  jw  is  substituted 
for  s  and  the  axes  represent  the  real  (Re)  and  imaginary 
(Im)  parts  o-f  G(jw)  -for  various  -frequencies. 
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A  Nyquist  Plot  is  given  in  Figure  IV. 2  and  it  can  be 
used  to  design  -feedback  gains  that  ensure  a  robust  system. 

The  measure  o-f  the  system  stability  are  Gain  and  Phase 
margin  obtained  from  the  Figure,  as  indicated. 


perturbed 
system 


/  +  Gfj(^J 


Re 
G(jU)) 


Figure  1^.2   Nyquist  plot  -  Stable  System 

The  Phase  Margin  (PM>  is  the  additional  phase  lag  at  the 
gain  crossover  -frequency  required  to  bring  the  system  to  the 
limit  o-f  instability  and  the  Gain  Margin  (GM)  is  the 
reciprocal   o-f   the  magnitude  |G(jm)  I  at  the  -frequency  where 
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the  phase  angle  is  -ISO  .  The  gain  crossover  -frequency  is 
the  -frequency  at  which  |G(jw)  I  is  unity. 

The  Phase  and  Gain  margins  are  a  measure  o-f  how  close  a 
polar  plot  is  to  the  -1  ■*■    jO  point. 

Considering  the  system  subject  to  additive  perturbation, 
as  we  have  in  the  diagram  o-f  Figure  IV.  3. 


U 


^O^ 


A  G 


Figure  IV. 3   Additive  Perturbation 

In  order  to  have  the  system  stable,  despite  the 
perturbation,  we  need 

lAG(jw)  I  <  II  +  G(jw)  I         (IV. 1) 

As  we  can  see  in  Figure  IV. 2,  this  condition  ensures  a 
stable  system. 

This  idea  can  be  extended  to  MIliO  problems  through  the 
use  o-f  matrix  norms  and  o-f  applying  the  Multi variable 
Nyquist  Theorem. 
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C.   MULTI-INPUT  MULTI-OUTPUT  SYSTEM 

The  generalization  o-f  the  SISO  theory  discussed  in  the 
previous  section  has  been  made  for  the  MIMO  problem. 

Let's  consider  our  system  as  represented  in  chapter  III 
and  apply  the  transformation  in  order  to  have  the  open  loop 
transfer  function  betMeen  input  and  output. 

X  =  A  X  +  B  u  (IV. 2) 

Applying   Laplace   Transform   to  the  equation  <IV.2),  we 


1 


have 


sX  =  AX  +  BU 


(IV. 3) 


wh 


ere   X  and  U  represent  the  Laplace  Transform  of  dx  and 


u(t),  respectively. Taking  the  value  of  X  in  equation  (IV. 1) 

,-1, 


As 

We  have 


X  =  Csi  -  a:  B  U 

Y  =  C  X 

Y  =  C  Csl  -  A:~^B 
U 


(IV. 4) 
(IV. 5) 

(IV. 4) 


This   transfer   function   corresponds   to  G(s)  in  Figure 


IV. 4. 


Figure  IV. 4   Return  Difference 
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From  Figure  IV. 4  it  is  possible  to  define  the  "Return 
Difference  Matrix"  (Kwakernaak  and  Sivan,  CRef.6D)  that 
corresponds  to  opening  the  loop  at  a  given  point  and 
connecting  an  external  input  variable  at  this  point;  for 
example,  at  point  1.  Assuming  the  input  u  as  zero,  we  have: 

Y^(s)  =  -  G(s)  H(s)  V  (s)        (IV. 5) 

where   V   (s>   corresponds  to  the  new  input  variable  and 

Y^(s)  the  "returned  variable". 

The  difference  between  Y  (s)  and  V(s)  is 

V(s)  -  Y^(s)  =  C  I  +  G(s)H(s)D  V(s)    (IV. 6> 

The  matrix  I  +  G(s)H(s)  is  defined  as  the  "Returned 
Difference  Matrix",  with  G(s)  equal  to 

G(s)  =  CCsI  -  A:~^B  (IV. 7) 

As  explained  by  Lehtomaki ,Sandel 1  and  Athans  in  CRef.7], 

the    Multi variable   Nyquist   Theorem   is  derived   from  the 

relationship 

detCI  +  G(s)H(S)D  =  Icl (s)       (IV. 9) 

iol  (s) 

where  $cl (s)  =  det (si  -  A)  corresponds  to  the 
characteristic  polynomial  of  the  open  loop  transfer  function 
and  lol (s)  =  det (si  -  A  +  BF)  corresponds  to  the 
characteristic  polynomial  of  the  closed  loop  transfer 
function  from  the  system  represented  in  Figure  IV. 4. 
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The  multi variable  Nyquist  theorem  which  requires  that  a 
closed  loop  stable  system  have  the  same  number  of 
counterclockwise  encirclements  of  the  origin  by  the  locus  of 
the  determinant  of  I  +  G(s)H(s)  as  the  number  of  open  loop 
poles  that  are    unstable. 

If  I  +  G(s)H<s)  is  quasi -singular ,  a  small  change  in  G 
may  make  the  matrix  singular  ;  this  causes  the  detCI  + 
6(s)H(s)D  to  become  zero  and  the  Nyquist  encirclement  count 
to  change  indicating  an  unstable  system. 

Basically, in  order  to  analyse  the  robustness  of  the 
system  we  have  to  verify  how  close  the  return  difference 
matrix  is  to  being  singular  as  a  function  of  frequency  (w). 

The  natural  measurement  of  the  singularity  is  the 
minimum  "singular  value",  since  this  is  the  tightest  norm. 

The  singular  value  of  a  matrix  M  is  defined  as 

8.     =  ex.  (  m'^M>>^^*  (IV.  10) 

1      1 

l_l 
M    represents   the   conjugate   transpose   matrix   of   a 

generic   M  matrix  and  X.  any  eigenvalue  of  the  product  of  M 

times  M. 

If  the  minimum  singular  value  is  close  to  zero  the 
matrix  is  quasi -singular  and  the  system  is  not  robust. 

If  one  assumes 

ffd  +  6(jw)H(jw))  >  oCjj         (IV. 11) 
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then  it  can  be  shown  CRe-F.  73  that  the  gain  and  phase 
margins  -for  the  system  may  be  represented  by  equations 
(IV. 12)  and  (IV. 13) ,  i.  e.  , 

Gain  margin:         GM  =    1  (IV. 12) 


—  1  2 

Phase  margin:        PM  =  +  cos  LI    -    anQ    2  (IV.  13) 


Universal  gain  and  phase  margins  curves,  presented  by 
Gordon  CRe-f.  SH  and  repeated  in  Figure  IV.  5  allow  the 
designer  to  pick  a  singular  value  that  corresponds  to  a 
desired  stability  margin  -for  the  system  to  be  designed. 

From  singular  value  plots  the  we  will  be  able  to 
identi-fy  the  critical  frequencies  or  range  of  -frequencies 
where  the  singular  value  is  below  the  chosen  level  from  the 
Figure  IV. 5. 
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Figure  IV. 5   Universal  Gain  and  Phase  Singular  Value  Plat 
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D.   SINGULAR  VALUE  ANALYSIS 

Using  the  program  developed  by  Gordon  CRe-f.SD,  with 
small  changes  to  make  it  compatible  with  the  size  o-f  our 
system  and  to  input  the  matrices  more  easily,  we  calculate 
the  minimum  singular  value  o-f  the  return  difference  matrix, 
with  the  frequency  varying  from  O  to  50  rad/sec. 

The  program  utilizes  the  IMSL  subroutines  to  generate 
the  return  difference  matrix  as  well  as  to  made  the  singular 
value  decomposition,  given  the  minimum  singular  value  for 
different  frequencies. 

The  return  difference  matrix  depends  on  the  point  where 
we  open  the  loop  represented  in  Figure  IV. 5.  For  the 
explained  case,  considering  point  1,  we  have  the  "output" 
difference  matrix,  that  corresponds  to  I  +  G(s)H(s);  if  the 
loop  is  open  at  point  2  we  have  the  "input"  case  and  the 
return  difference  matrix  will  be  I  +  H(s)G(s). 

From  the  Figure  IV. 5  we  choose  as  reference,  for  our 
analysis,   the  minimum  singular  value  as  0.6  that  gives  us  a 

0 

phase  margin  of  about  35   to  characterize  a  "robust"  system. 

The  minimum  input  singular  value  is  plotted  in  Figure 
IV. 6. 

It  shows  that  the  system  is  robust,  except  for 
frequencies  between  0.6  and  5.0  rad/sec. 
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The  results  related  to  the  minimum  output  singular  value 
are  in  Figure  IV.  7  and  the  robustness  o-f  the  system  is  poor 
from  0  to  18  rad/sec. 

Those  results  indicate  that  the  behavior  of  the  system 
under  perturbations  should  present  problem  for  frequencies 
below  10  rad/sec. 

To  confirm  this  expectation,  we  have  plotted  the  BODE 
diagram  (magnitude)  for  the  closed  loop  system  considering 
the  angle-of -attack  and  altitude  with  respect  to  both  inputs 
(acceleration  and  desired  altitude). 

The  diagrams  are  on  Figures  IV. 8  to  IV. 11 ,  where  can  see 
the  attenuation  at  the  low  frequencies  indicating  that  we 
can  not  expect  good  response  at  those  frequencies  and  one 
could  further  expect  problems  with  robustness  as  is  evident 
from  the  singular  value  plots. 

On  the  next  chapter  we  try  to  improve  the  system 
designing  a  new  set  of  feedback  gains,  in  order  to  have  a 
robust  system. 
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SINGULAR  VALUE  ANALYSIS 
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Figure    IV. 6       Input    Minimum    Singular    Value 
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SINGULAR  VALUE  ANALYSIS 
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Figure  IV. 7   Output  Minimum  Singular  Value 
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Figure  IV. 9 
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Figure  IV. 11 
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V.  IMPROVING  THE  DESIGN 

Under  this  chapter  Me  try  to  improve  the  system  by 
calculating  new  feedback  gains  that  will  yield  a  more  robust 
system,  taking  into  consideration  the  results  indicated  in 
Figures  IV. 6  and  IV. 7. 

The  POPLAR  program,  developed  by  Gordon  CRef.SD  will  be 
used.  A  numerical  optimization  technique  is  applied  in  order 
to  increase  the  minimum  singular  value  o-f  the  considered 
return  difference  matrix,  therefore,  a  resultant  robustness 
of  the  design. 

The  development  in  this  chapter  will  be  preceded  by 
comments  on  numerical  optimization  and  a  description  of  the 
computer  program. 

A.   OPTIMIZATION 

The  optimization  was  accomplished  by  means  of  the 
Automated  Design  Synthesis  Program  (ADS)  developed  by 
Vanderplaats  CRef.9D. 

The  purpose  of  ADS  as  of  others  numerical  optimization 
routines  is  to  find  the  "best"  possible  solution  for  the 
problem,  starting  from  an  initial  set  of  variables  and 
updating  the  design  i teracti vely.  The  problems  can  develop 
in   convergence   of   the   method   and   in   the  computer  time 
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needed.  I-f  the  problem  has  multiple  solutions,  the 
optimization  does  not  always  lead  to  the  absolute  optimum. 

The  ADS  program  is  designed  as  a  black  box  optimizer 
which  allows  the  user  to  choose  combinations  o-f  one 
dimensional  search,  optimization  algorithm  and  strategy. 

ADS  is  used  as  a  subroutine  and  the  parameters  that 
correspond  to  the  di-f-ferent  applications  are  chosen  by  the 
user,  as  explained  in  CRe-f.9]. 

The  preferred  method  used  in  the  ADS  applications  is 
re-ferred  to  as  Sequential  Unconstrained  Minimization 
Techniques  and  can  be  considered  as  a  method  that  starts 
with  an  objective  function  and  the  constrains  combined  into 
an  augmented  objective  function  and  then  minimizing  this 
function  as  if  the  problem  was  unconstrained. 

ADS  employs  penalty  function  techniques  as  well  as  an 
Augmented  Lagrange  Multiplier. 

As  presented  in  CRef.BD,  the  iteration  between  the 
user's  program  and  the  ADS  routine  can  be  represented  by  the 
block  diagram  of  Figure  V.l;  where  the  "info"  parameter  is 
used  as  a  "flag"  for  a  dialogue  between  the  user  and  the  ADS 
routines. 

The  details  for  utilization  of  ADS  are  presented  by 
Vanderplats  in  CRef.SD. 
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DIMENSION  ARRAYS 


INITIALIZE  PARAMETERS 


Yes 


INFO  =  O 


CALL  ADS 


EVALUATE 
OBJECTIVE 
&  CONSTRAINTS 


No 


/\r 

HU3 

CONTROL 

" 

STRATEGY 

■  " 

OPTIMIZER 

'• 

ONE 
DIMENSIONAL 
SEARCH 

Yes 


EXIT 
SOLUTION 
COMPLETE 


EVALUATE 
REQUIRED 
GRADIENTS 


Figure  V.  1   Organization  o-f  ADS  Program 
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B.   POPLAR  PROGRAM 

The  Gordon's  program  has  the  capability  of  calculating 
the  -feedback  gains  that  improve  the  response  of  the  system 
either  by  pole  placement  and/or  by  increasing  the  minimum 
singular  value  of  the  return  difference  matrix. 

The  main  objective  is  to  improve  the  robustness  by 
increasing  the  minimum  singular  value  above  a  desired  level 
chosen  from  the  universal  curve  reproduced  in  Figure  IV. 5. 

"Optimum"   values   of   the   feedback   are   calculated  by 

minimization   of   an  Objective  function  whose  pole  placement 

part  is 

1 

""^^    =X^  >^R-    ^R.^'  "  '    S  -  ^I.>* 
^^  ^   di    1  di    1 

1  =  1 

Mhere 

X_   -  real  part  of  the  desired  eigenvalue  i; 
R  .  . 
dx 

X    -  real  part  of  the  computed  eigenvalue; 
R . 

1 

X    -  imaginary  part  of  the  desired  eigenvalue  i; 

di 
X    -  imaginary  part  of  the  computed  eigenvalue  i; 

i 

and,   for   the   minimum   singular  value  optimization,  we 
have 

OBJ  =  J>  <:  max  CO,  «.  -ff.^JW,p)  3  3 


where  a         indicates   the  desired  minimum  singular  value 
and  ^  is  the  minimum  singular  value  at  a  certain  frequency. 
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The  optimization  procedure  changes  the  -feedback  gains 
until  the  minimum  singular  value  is  raised  to  the  desired 
level . 

The  pole  placement  and  robustness  program  calculates  the 
return  difference  matrix  to  be  considered,  working  in  the 
complex  space  (as  the  analysis  is  conducted  in  the  -frequency 
domain)  and  the  objective  -functions,  calling  the  ADS  routine 
to  make  the  optimization. 

Input's  -for  the  program  are  the  matrices  o-f  the  state 
representation  o-f  our  missile  (  A  ,  B  ,  C  and  F), where  the 
observation  matrix  (C)  was  specified  as  Identity  because  we 
are  assuming  all  the  states  are  observable  as  well  as 
control lable. 

Other  inputs  are  the  desired  minimum  singular  value, 
desired  pole  locations  and  frequency  interval  to  be 
considered. 

In   our   application  of  POPLAR  the  initial  values  of  the 

feedback   gains  are      those   of   the  original  system  and  the 

parameters  for  application  of  the  ADS  program  are: 

Strategy  -  Augmented  Lagrange  multiplier; 

Optimizer  -  BFGS  variable  metric  method  for  unconstrained 
minimization; 

One— Dimensional  Search  -  Polynomial  interpolation. 

Two   situations   will  be  considered;  the  INPUT  -  minimum 

singular  value,  where  the  return  difference  matrix  is  I  +  HG 

and   OUTPUT-minimum  singular  value  with  I  +  GH  as  the  return 

difference  matrix. 
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C.   INPUT  -  MINIMUM  SINGULAR  VALUE 

The  program  was  first  run  to  improve  the  minimum 
singular  value  -for  the  input  situation  within  the  same 
interval  of  frequencies  used  in  Chapter  IV. 

Due  to  the  large  CPU  time  involved,  the  program  was  run 
for  the  critical  band  of  frequencies  and  the  result 
evaluated  for  the  complete  interval. 

The   minimum   singular  value  taken  as  reference  from  the 

0 

Universal  curve  was  0.6  that  implies  a  phase  margin  of  35  . 

We  have  an  objective  function  taking  into  consideration 
the  improvement  of  the  minimum  singular  value  as  well  as  a 
part  corresponding  to  the  pole  placement.  The  main  goal  is 
to  improve  the  singular  value  but  with  a  solution  such  that 
the  poles  are  in  positions  where  the  time  response  satisfies 
the  requirements. 

The  weight  of  the  pole  placement  part  was  considered  as 
107.    of  the  singular  value  part. 

The  program  was  started  considering  only  part  of  the 
feedback  matrix  as  free  parameters,  increasing  this  number 
until  a  reasonable  solution  was  found. 

The  best  solution  was  found  keeping  the  original 
feedback  gains  and  calculating  gains  to  feedback  all  the 
other  states  to  generate  the  commanded  acceleration. 

The  computer  output  of  the  best  solution  with  the 
corresponding  inputs  is  given  in  Appendix  B. 
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The   -feedback   gains   that  improve  the  robustness  o-f  the 
system  are  the  -following 


F    = 


0.04    -1.15    0.03    -1.7    -9.07    0.32    7.07    0.049    0.986    1.4 
000  0000  O  00 


The  plot  o-f  the  minimum  singular  value  -for  dif-ferent 
-frequencies  is  given  in  the  Figure  V.2,  along  with  the 
values  -for  the  original  system. 

As  we  can  see  the  minimum  singular  value  was  increased 
in  the  critical  values  of  frequency  with  the  minimum 
changing   from  0.40  to  0.63  that  indicates  an  improvement  in 

0  0 

the  phase  margin  from  20  to  35  .  Also  the  bandwidth  where 
the  system  was  less  robust  was  changed  from  5.0  rad/sec  to 
approximately  1.0  rad/sec. 

The  robustness  of  the  control  system  for  all  the 
considered  frequencies  is  better,  arriving  close  to  0.9  with 
a  considerable  improvement  at  low  frequencies. 

Using  the  software  Controls,  the  time  response  was 
plotted  in  Figure  V.3. 

The  maximum  and  minimum  values  of  the  angle-of-attack 
are  smaller  than  those  of  the  original  system;  the  maximum 
acceleration  and  velocity  are  slightly  lower  and  an 
undershoot  of  about  lOX  appears  on  the  altitude,  that  could 
be  a  problem  in  some  flight  conditions. 
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Me  now  turn  to  a  consideration  of  the  output  return 
difference  matrix  although  this  is  less  important  from  a 
practical  viewpoint. 
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Figure    V.2   -    Improvement    in    the    INPUT-Minimum   Singular    value 
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Figure  V.3  -  Time  Response  -for  the  Improved  System  (INPUT) 
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D.   OUTPUT  -  MINIMUM  SINGULAR  VALUE 

In  the  Figure  V.4  we  have  plotted  the  output  minimum 
singular  values  that  correspond  to  the  run  where  the  input 
singular  values  were  improved,  as  well  as  the  original 
values. 

As  a  considerable  degradation  o-f  the  output  singular 
values  was  evident,  we  try,  in  the  last  part  o-f  this  work, 
to  arrive  to  a  situation  where  the  system  could  be 
reasonable  robust  in  both  cases,  increasing  the  output 
singular  value  even  with  some  reduction  at  the  input 
situation. 

After  several  runs  of  the  program,  using  different 
combinations  of  optimization  techniques  with  different 
starting  points,  it  was  verified  that  due  to  the 
characteristics  of  the  plant,  that  increases  in  the  input 
singular  values  resulted  in  decreases  in  the  output 
situation. 

In  order  to  obtain  some  improvement,  an  effort  was  made 
to  modify  the  system. 

The  first  step  was  to  look  at  the  controllability  matrix 
of  our  model ,  represented  by 

1.855632D-01    4.948168D+02 

2.725664D+13   -2. 690449D+13 

2.725664D+13   -2. 690449D+13 

2.8611210+00   -2.822009D+00 

2.201811D+01    2.171710D+01 
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-6.7S1791D+00    6.6891110+00 

-2 . 286975D+00    2 . 255563D+00 

9 . 855433D+00   -9 . 7 1 96 1 4D+00 

- 1 . 508382D+0 1    1 . 487685D+0 1 

4.905487D-01   -4. 264398D-01 

As  we  can  veri-fy,  all  the  elements  are  non  -ero,  this 
means  that  the  poles  can  placed  at  any  desired  position;  but 
a  close  look  at  the  2nd  and  3th  rows  indicates  very  high 
order  numbers  compared  to  the  others.  This  situation 
generates  some  numerical  difficulties  in  pole  placement. 

Using  engineering  judgement,  we  have  to  de-Fine  what  kind 
o-f  changes  have  to  be  made  -for  achieving  a  more  robust 
system. 

First  we  try  to  get  a  balanced  A-matrix  by  changing 
units,  i.  e.  ,  the  angles  will  be  in  radians  instead  o-f 
degrees  and  the  pitch  angular  rate  in  radians  per  second 
instead  o-f  degrees  per  second.  Despite  some  reduction  on  the 
numbers  was  not  possible  to  arrive  to  a  satis-f  actory 
solution. 

More  positive  effect  was  obtained  by  adding  the  effect 
of  the  angle-of-attack  in  the  controller  (see  Figure  II. 3). 

As  the  original  system  has  the  output  singular  values 
higher  than  0.5  for  frequencies  above  10  rad/sec,  we  have 
kept  the  original  feedback  gains  and  calculated  gains  to 
feedback  only  the  angle-of-attack,  pitch  rate  and  commanded 
actuator. 
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The  gains  are  the  following: 


F    = 


0  0  0  0.037       0.027       O       -0.02       0.049    0.986    1.404 

000  O  00  O  O  O  O 


The  result  is  presented  in  Figure  V.5,  compared  with  the 
original  singular  values  and  those  from  the  improvement  in 
the  input  case. 

The  output  singular  values,  compared  with  the  original 
was  slightly  improved  for  frequencies  below  5.0  rad/sec. 

Resultant  input  singular  values  Are  plotted  in  Figure 
V.6,  and  they  are  lower  than  those  from  the  improvement  in 
the  input  situation  but  a  little  higher  than  the  original. 

The  time  response  were  practically  unchanged  with 
respect  to  the  original  model  and  it  is  presented  in  Figure 
V.7. 

Further  analysis  is  needed  to  obtain  a  solution  with 
high  output  singular  values. 
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OUTPUT  -  MINIMUM  SING.  VALUE 
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Figure  V.4   Output  Singular  Value  Resulted  from  the 
Improvement  in  the  Input  Singular  Value 
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OUTPUT  -  MINIMUM  SING.  VALUE 
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Figure  V.5   Output  Singular  Value 


103 


INPUT  -  MINIMUM  SING.  VALUE 
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Figure  V.6   Final  Result  in  the  Input  Singular  Value 
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Figure  V.7   Final  Time  Response 
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VI.  CONCLUSIONS 

The  sensitivity  analysis  based  in  the  change  in 
eigenvalues  represents  an  important  role  in  the  design. 

As  was  demonstrated  -for  the  studied  missile  its  behavior 
corresponds  to  a  system  insensitive  to  the  variation  o-f  the 
considered  aerodynamic  parameters. 

The  analysis  gives  a  good  indication  o-f  Mhat  parameters 
have  to  be  precisely  determined. 

The  robustness  analysis  also  demonstrated  how  useful  the 
singular  value  analysis  is  as  an  auxiliary  tool  for  the 
designer. 

Singular  value  analysis  indicates  that  the  system  is 
robust  to  input  perturbations  but  is  deficient  in  robustness 
to  output  perturbations. 

The  physical  nature  of  the  problem  indicates  that  the 
prime  concern  should  be  with  respect  to  the  input 
perturbations;  the  output  perturbations  are  of  minor  concern 
but  should  be  kept  in  mind  if  unusual  conditions  should  be 
encountered  by  the  missile. 

On  the  MIMO  design  a  commomly  used  method  is  the  Linear 
Quadratic  analysis  where  the  performance  levels  are  reached 
by  adjusting  weighting  terms  in  the  "cost"  function,  but  the 
results  of  that  method,  for  non-diagonal  R  matrices,  do  not 
necessarily    imply    a   robust   system.   Using   the   return 
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di-f-ference  matrix  Me  can  improve  the  robustness  o-f  the 
design. 

One  advantage  o-f  this  procedure  is  that  it  permits  a 
high  level  o-f  interaction  between  the  system  and  the 
designer. 

Further  problems  of  a  computer  nature  (i.e.,  large  CPU 
time)  are  encountered  with  a  high  number  o-f  states.  Analysis 
must  be  con-fined  to  a  small  frequency  range  of  low  singular 
values  and  the  complete  system  return  matrix  singular  value 
are  calculated  with  the  determined  feedback  gains. 

One  aspect  that  should  be  considered  in  future 
development  is  to  improve  the  cpu  time  used  when  handle  a 
system  with  a  high  number  of  states. 
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APPENDIX  A 
MISSILE  DATA 
Under   this  Appendix  the  missile  sizing,  mass  properties 
and   aerodynamic   parameters  are  presented  as  given  by  Arrow 
CRe-f.l],  including  the  figures. 

The  missile  is  1/6  scale  o-f  the  actual  circular  missile 
configuration  and  is  reproduced  in  Figure  A.l.  It  is  tail 
controlled   using   four   identical   control  surfaces  located 

0 

with  +30   dihedral. 


A.   GEOMETRY  AND  MASS  PROPERTIES 

In  table  VI  we  have  the  size  and  mass  properties,  with 
the  respactive  values,  used  for  development  of  the  state 
equations. 

Only  the  uncoupled  pitch  channel  was  considered, 
assuming  no  roll  movement. 


Table  VI 
MISSILE  GEOMETRY  AND  MASS  PROPERTIES 


Weight 
Length 
Diameter 

Reference  Area  for  Coeff. 

Moment  of  Inertia  about  y 


M 
1 
d 


S 
I 


yy 


2525  lbs 

168   in 

24   in 

2 
IT  ft 

804  slug-ft^ 
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Figure  A.  1   Circular  Configuration  o-f  the  Missile 
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B.   AUTQPILOT/AIRFRAME 

Figure  A. 2  shoMS  the  block  diagram  of  the  autopilot  used 
to  develop  the  model  o-f  the  pitch  control  channel  of  the 
missi le. 

The  state  equations  that  correspond  to  the  autopilot  are 
the  following: 

x^=    -150.0    X       +    150.0    D 

X    =    1 .  353    X       +    X       -    1 .  353    i) 

2  13  'zc 

X    =    -6.572    X       -    5.0    X       +    6.572    i) 

3  13  'zc 

X    =    -44.3316    X       -    59. 11    X 

H  5  fe 

x_=    X       -    0.14«2    x^    +    0.0395    x^ 
x^=    -188.4    X       +    188.4    x_ 

O  O  / 

X  =  -0.4608  X   -  2.231  x   -  0.3406  x   +  2.231  x 

7  12  3  4 

-  15.095  x^  -  20.13  X   -  0.1430  x   +  0.4608  d 

S  fr  7  ZC 

The    input    for    the    autopilot    is   the   commanded 

acceleration i) 

zc 

The  values  related  to  the  aerodynamic  are  the  following: 


c 

mo: 

^ 

-    0.06 

C 
msp 

=    - 

-    0.08 

So: 

= 

0.  15 

^N&p 

= 

0.04 

110 


—  2 

Dynamic  Pressure   -   q   1650     lb/ft 

Velocity  -   V   3825.46  ft/sec 
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APPENDIX  B 
POPLAR  PROGRAM 
The      computer       output    where   the    improvement    in    the    input 
singular    value    was    found    is    listed    below. 
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